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<^ , Abstract. We apply elliptical isothermal mass models to reproduce the 

| point source properties, i.e. image positions, flux density ratios and time 

■ delay ratios, of the quadruple lens B1608+656. A wide set of suitable 

<D ' solutions is found, showing that models that only fit the properties of 

point sources are under-constrained and can lead to a large uncertainty 
in the prediction of H D . We present two examples of models predicting 
H = 100km s^Mpc -1 (x 2 = 4) and H Q = 69kms- 1 Mpc~ 1 (x 2 = 24). 
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, 1. Introduction 

C^l | Relative positions, flux ratios and time delays of the 4 images in B1608+656 have 

been presented here by Fassnacht (1999) and references therein (c/ Table ||). 
Koopmans & Fassnacht (1999) have concluded H Q = 59lgkms _1 Mpc _1 within 
the context of a family of parametrized, isothermal models. Here, we investigate 
whether a larger set of models allows a wider range of Hubble constants. 
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£3 . 2. Elliptical isothermal models 

Following Blandford & Kundic (1997), we adopt a scaled lensing potential ifj 
composed of two elliptical contributions to describe the lensing galaxies Gl and 
^ I G2 plus external shear 7: 

V>=^ bi { S f + rf [1 - a cos(2(pi - fa))}}^ + r 2 7 cos(2( V 9 1 -^ 7 )) (1) 
i=i 

Here (rj, ^) are polar coordinates with origin at the center of each galaxy, s 
measures the core radius, e and <f> the ellipticity and position angle of the major 
axis. At large radius the mass distribution is isothermal, going as X oc r . The 
lenses will be fixed at x Gl = (0.446, -1.063)" and x G2 = (-0.276, -0.937)", the 
centroids in H band, which are less affected by reddening (Blandford, Surpi &: 
Kundic 1999). 

We minimize a x 2 function. The best fit achieved, hereafter Model A, has 
X 2 = 4.0 and yields H = 100km s _1 Mpc _1 . Models with lower values of H a 
can also be built fixing H Q and fitting the 3 time delays instead of the time 
delay ratios. As an example we present the results of Model B having H = 
69kms -1 Mpc~ 1 and x 2 = 24.7. The parameters and predictions of Model A and 
B are displayed in Tables p] and || respectively. They represent reasonable mass 
distributions given, especially, our ignorance of the dark matter distribution 
(Figure 1). 
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Figure 1. Surface mass density in Models A and B 
3. Discussion 

A variety of parametrized models can reproduce the point source properties of 
B1608+656. This precludes an accurate determination of H Q . To break the de- 
generacy, extra constraints, associated with the extended emission of the source, 
have to be incorporated. It is also helpful to specify the distribution of dark mat- 
ter on larger scale than the image distribution. A similar conclusion has been 
drawn by Williams & Saha (1999) using pixellated models. 



Table 1. 


Model parameters. 










Model A (Ho = 100km s^Mpc- 1 ) 


Model B (Ho = 


69kms~ 1 Mpc~ 1 ) 


Parameters 


Gil ii'2 


Gl 


ii'2 


s(" 


) 


0.10 0.10 


0.00 


0.05 


b 


0.9072 0.2453 


0.7797 


0.3429 


e 




0.3269 0.6405 


0.1570 


0.3149 


<t>(° 


) 


163.45 154.93 


172.62 


160.86 




a 


0.0876 -10.92 


0.0473 


12.47 


Table 2. 


Comparison between observat 


ions and model predictions. 


Properties 




Observation 


Model A 


Model B 


*a(") 


( 0.0000, 0.0000) ± (0.0023,0.0023) 


( 0.0000, 0.0000) 


( 0.0000, 0.0000) 


xb(") 


(-0.7380,-1.9612) ± (0.0043,0.0046) 


(-0.7382,-1.9613) 


(-0.7365,-1.9518) 


xc(") 


(-0.7446,-0.4537) ± (0.0045,0.0049) 


(-0.7443,-0.4544) 


(-0.7422,-0.4575) 


&(") 


( 1.1284,-1.2565) ± (0.0107,0.0124) 


( 1.1271,-1.2582) 


( 1.1269,-1.2207) 


Fa/Fb 




2.042 ± 0.124 


1.917 


1.901 


Fc/Fb 




1.037 ± 0.083 


1.092 


1.131 


F d /Fb 




0.350 ± 0.055 


0.428 


0.504 


Tab (d) 




26.0 ± 5.0 


28.4 


27.6 


Tcb (d) 




33.0 ± 5.0 


32.0 


31.7 


T DB (d) 




73.0 ± 5.0 


68.4 


71.3 


x 2 




0.0 


4.0 


24.7 
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